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We study the running of the QCD coupling with the momentum squared (Q 2 ) and the tempera- 
ture scales in the high temperature limit (T > T c ), using a mass dependent renormalization scheme 
to build the Renormalization Group Equations. The approach used guaranty gauge invariance, 
through the use of the Hard Thermal Loop approximation, and independence of the vertex chosen 
to renormalize the coupling. In general, the dependence of the coupling with the temperature is not 
logarithmical, although in the region Q 2 ~ T 2 the logarithm approximation is reasonable. Finally, 
as known from Debye screening, color charge is screened in the coupling. The number of flavors, 
however, is anti-screened. 
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£Nl ' One of the important questions in perturbative QCD at finite temperature T is what is the temperature dependence 
of the strong coupling constant, a s (T). According to an early study by Collins and Perry [1], where arguments from 
• the Renormalization Group Equations (RGE) were used, temperature replaces the momentum as the running scale, 
| implying that a s (T) decreases as a logarithm as T increases. This conclusion is based on general considerations, 
where no specific calculations of the vertices and self-energies were done. Later, more rigorous treatments based on 
| the RGE were developed [2,3], with applications to perturbative QCD [2,4]. However, it was soon realized that the 
resulting coupling was strongly dependent on both the gauge [5-7] and vertex chosen to renormalize it [4,8,9]. In 
t-H ' any case, some of the calculations did not result in the logarithmical dependence of the coupling in the temperature. 
Recently, in the context of a semiclassical background field method, it has also been calculated [10] that a s has an 
inverse logarithmical dependence in the temperature, where reference to the earlier calculations based on the RGE 
was said to substantiate the result. Moreover, applications to the perturbative calculations of energy density, pressure, 
i etc, using the logarithmical dependence are frequently used in the literature [11-13]. 

In the early 1990s the problem of gauge invariance of two, three and four points Green functions was solved with the 
introduction of the Hard Thermal Loops (HTL) [14]. In particular, it was shown that the dominant, gauge invariant, 
contributions to the N-point functions have a T 2 dependence, and that they satisfy abelian type Ward identities 
[14,15]. From the gluon self-energy in the HTL, it is possible to build an effective coupling at finite temperature, which 
goes with the inverse of the Debye mass and thus does not have a logarithmic dependence on the temperature [16]. 
However, there has been no extensive calculation of the RGE for the ultra-violet regime of the running coupling using 
the dominant, gauge invariant, results from the HTL. Rigorously, there is no extra ultra-violet divergence induced by 
thermal effects, meaning that the renormalization of QCD at T — is sufficient to make QCD at T / finite. Still, we 
• • can redefine the theory by a finite renormalization, relating parameters like mass and coupling constants at different 
temperatures with the aid of temperature dependent renormalization constants, with appropriate renormalization 
conditions [3]. The purpose of this letter is to calculate the running of the QCD coupling with the momentum scale 
and with the temperature using the vertices and self-energies calculated from the HTL. Although it is an immediate 
application, we think that given the discrepancies found in the literature, this is an useful exercise with an impact in 
all the perturbative QCD calculations at finite temperature. 

We will treat the mass scale introduced in the system by the temperature in the same way quark masses are 
handled in the T = theory [17], in a mass dependent renormalization scheme. In this case, the general form for the 
renormalized n-points Green function will be: 

T n =Y n {k,g 2 ,k 2 /y 2 1 T 2 /^ 2 ), (1) 

where fi is the renormalization scale, k is the external momentum and g the renormalized coupling at the scale /i. 
Using the standard approach of the RGE, i.e., taking the derivative of the unrenormalized Green function with respect 
to fi, one gets: 
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where 7„ ^g 2 , is the anomalous dimension of the n fields appearing in T n and 



r„ = o, (2) 
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the beta function and the temperature anomalous dimension, respectively. The temperature dependence in the beta 
function and anomalous dimensions indicates that we are working in a mass dependent renormalization scheme. The 
relation between the bare and renormalized coupling is given by: 



a s (» 2 ,T 2 ) = Z- 1 ( fJ ?,T 2 )a , 



(4) 



where Z a is the renormalization constant of the coupling, ao = 5o/4tt is the bare coupling, and a s — ,g 2 /47r is the 
coupling at the subtraction point /i and temperature T. 

If we want a meaningful result for the running coupling, any vertex should be equally good to renormalize it. If we 
use the triple gluon vertex, then 



z-v,t7a30 



Zj^ 2 ,T 2 /A 2 T ) 



(5) 



where Z% renormalizes the gluon field while Z\ rcnormalizcs the triple gluon vertex. We start with the renormalization 
of the gluon field. To include temperature in the renormalization constants and simultaneously preserve the Slavnov- 
Taylor identities, we will work in the framework of the HTL resummation program, where the dominant terms are 
known to be gauge invariant. A direct computation of the transverse part of the gluon self-energy at one loop gives 
[18,19]: 
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where k is the four momentum of the external gluon. Only the dominant term in the high temperature expansion was 
written because of its gauge invariance. One can see that Eq. (6) vanishes for ko << \k\, and that it is reduced to 
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for ko » \k\, and At is a mass scale (A T = \k\ko). 

We will be working in a Momentum Subtraction Scheme (MOM), which is suitable for a mass dependent renormal- 
ization scheme as in the present case, with the temperature replacing the mass. Although a MOM scheme usually 
breaks the Slavnov- Taylor identities, the deviation among the couplings defined through different vertices is vanish- 
ingly small [20] . As a renormalization condition, we impose that the thermal part of the renormalized gluon self-energy, 
calculated at the subtraction point /x 2 , be given by: 



n 



(i) 

T,R 



n f \ 2tt 2 T 2 



,2 ' 



2 / 3 /r 

Including the T — part, we then calculate the renormalization constant Z 3 : 



z, = i + z^ = i + ^( Nc + n -l) 
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where a is the gauge parameter. We now need the one loop correction to the triple gluon vertex, which is a known 
quantity [16]: it has the same functional form of Eq. (7). Imposing again that the temperature dependent part of the 
renormalized vertex be equal to the bare one, with the replacement of At — ► we calculate the triple gluon vertex 
renormalization constant to be: 



Zl = l + Z^ = l + ^( 
Using Eqs. (9) and (10) in (5), we get: 
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Notice that the temperature dependence, unlike the momentum scale, is not logarithmical. Before we proceed to 
calculate the running of the coupling, we make a few remarks on the calculation of Z^ 1 using other vertices. For 
instance, if we use the ghost-gluon vertex, then: 



Z^Zz 



(12) 



where Z 3 is the ghost field renormalization constant and Z\ is the gluon-ghost vertex renormalization constant. It 
happens that in the HTL approximation the external momenta in the numerators of the loop integrals are disregarded, 
as the main contribution comes from the internal momenta of order T, with T taken to be large. Compared to the gluon 
self-energy, the ghost self-energy does not have enough powers of internal momenta in the numerator to produce the 
leading, T 2 , behavior in the temperature. The same happens to the ghost vertex, something that can be immediately 
inferred from the abclian Ward identities relating the ghost vertex to its self-energy. Hence, the thermal parts of Z% 
and Z\ will be sub-leading to Z?, in the temperature. In such case, the thermal part of Z a will be dominated by Z3 
and we recover Eq. (11). Finally, we could have used the quark-gluon vertex to calculate Z a : 
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where Z 2 is the fermion field renormalization constant while Zip is the quark-gluon vertex renormalization constant. 
The HTL calculation for the dominant part of the quark self-energy and the quark-gluon vertex at one loop are 
known [16]. Their T 2 dependence are the same. Imposing renormalization conditions of the type of Eq. (8) for the 
calculation of Z 2 and Zip, it follows that their T 2 dependence are the same, cancelling in the ratio of Eq. (13). As 
before, Z 3 will dominate the temperature dependence of Z a , and we once again recover Eq. (11). 

The calculation of the (3 function (including the T = and the T / parts), for a fixed temperature but an 
arbitrary renormalization point is now straightforward. Using Eqs. (3), (4), and (11), we have: 
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to order a 2 . This is the only RGE for a s because there is only one renormalization scale for both the T = and 
T ^ parts. The solution of Eq. (14) is: 
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It is helpful to rewrite Eq. (15) in the same format of the T = theory. To this end, we define effective, scale 
dependent, numbers of colors and flavors: 
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With the help of Eqs. (16) and (17), the expression for the running coupling is written as: 

a s (QlT 2 ) 
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Equations (16) and (17) tell us that for fixed Qo an d T 2 , with Q 2 — > 0, the effective number of color decreases, 
N*ff < N c , while the effective number of flavors increases: color charge is screened, and the number of flavors are 



anti-screened. The same argument applies when keeping Q 2 and Q 2 fixed, while taking T 2 going to infinity. To 
quantify these assertions, we show in Fig. 1 the Q 2 dependence of rij and N*ff for T 2 = 1 GeV 2 , Qq = m|, and 
N c = rif = 3. As expected, N^l decreases while n e P increases as Q 2 decreases from Ql to T 2 . 
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FIG. 1. The effective number of colors and flavors at T — 1 GeV as a function of Q 2 . For large values of Q 2 , both tend to 
their values in the T = GeV theory. 

In Fig. 2 we show the behavior of a s for 3 values of the temperature. For a s (Ql,T 2 ), we use the experimental 
value measured at mz ~ 91 GeV and at zero temperature [22]. That is, we assume that at such high values of the 
virtuality of the probe, temperatures of the order of 1 GeV are not relevant. This seems to be the case as, in this 
region, N^ fi — > N c and n e ^ — ► «/. At T = GeV we have, as usual, the coupling growing rapidly for Q 2 < 10 GeV 2 . 
However, at T = 0.5 and 1 GeV, a s (Q 2 , T 2 ) starts to change its behavior in the region around Q 2 = 20 GeV 2 . Instead 
of the rapid growth observed at the T = case, for finite T there is first an almost Q 2 independence of the coupling, 
and then it decreases with Q 2 , showing the color screening, and flavor anti-screening, in action. In fact, we observe 
that a s changes its qualitative behavior (from a divergence to finite values), for small values of Q 2 , in the region 
around T ~ 0.2 GeV, although any conclusion based on a perturbative RGE analysis for such low values of T and 
Q 2 should be taken with extreme caution. In any case, in general the running coupling does not have a logarithm 
dependence on the temperature. But if we consider the particular region where Q 2 ~ T 2 , taking Qq >> T 2 , we get 
from Eq. (16-18) that: 

a s (Q 2 ~ T 2 ,T 2 ) - a s (T 2 ) x ^ (19) 

{llN c -2nf)ln[ 7 ^-) 
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FIG. 2. The strong coupling constant as a function of Q 2 calculated for 3 differerent values of the temperature. 
In summary, we have presented a RGE calculation of the running coupling at finite temperature QCD, introducing 
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the temperature scale as the mass scale is introduced in the T = theory, in a mass dependent renormalization 
scheme. We restricted the calculation to the dominant, gauge invariant, HTL terms. As a result, the coupling is 
shown to be independent of the vertex used to renormalize it. As is well known from Debye screening, the number of 
colors are screened by the temperature. The number of flavors, however, is anti-screened. Also, although in general 
the temperature dependence of the coupling is not logarithmic, in some especial cases, where Q 2 ~ T 2 , the logarithmic 
dependence can be a good approximation. 
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